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Abstract 

We propose a new mechanism to generate density perturbations in inflationary mod- 
els. Spatial fluctuations in the decay rate of the inflaton field to ordinary matter 
lead to fluctuations in the reheating temperature. We argue that in most realistic 
models of inflation the coupling of the inflaton to normal matter is determined by 
the vacuum expectation values of fields in the theory. If those fields are light dur- 
ing inflation (this is a generic situation in the minimal models of super symmetric 
inflation) they will fluctuate leading to density perturbations through the proposed 
mechanism. We show that these fluctuations could easily dominate over the ones 
generated through the standard mechanism. The new scenario has several conse- 
quences for inflation model building and observations. The proposed mechanism 
allows to generate the observed level of density perturbations with a much lower 
scale of inflation and thus generically predicts a smaller level of gravitational waves. 
The relation between the slope of the spectrum of the produced density perturba- 
tions and the potential of the inflaton field is different from the standard relations 
obtained in the context of slow roll inflation. Because the field responsible for the 
fluctuations is not the inflaton, it can have significantly larger self couplings and 
thus density perturbations could be non-Gaussian. The non-Gaussianity can be 
large enough to be detectable by CMB and Large Scale Structure observations. 



1 Basic Mechanism 



In the standard picture [T], the observed density perturbations are produced as 
follows. As the inflaton field <fi rolls down its potential it is effectively massless, 
so it fluctuates up and down. These fluctuations are different in different regions, 
resulting effectively in inflation lasting different amounts of time in different places. 
At the end of inflation, the inflaton oscillates about the minimum of its potential 
and decays, reheating the universe. As a result of the fluctuations each region of 
the universe goes through the same history but at slightly different times leading to 
adiabatic density perturbations. 

Our idea is different. To reheat the universe the inflaton has to couple to the 
ordinary particles through a coupling schematically given by, 



where q stands for the ordinary particles (e.g. quarks, leptons and their scalar 
superpartners), and A is the coupling strength, which is assumed to be constant in the 
standard picture. For simplicity we assume that the density perturbations produced 
at the inflationary stage are negligible. So inflation ends at the same time everywhere 
in the universe. The density fluctuations are created during the reheating process, 
because of fluctuations in the coupling "constant" A. Since A controls the efficiency 
by which the energy stored in the 0-condensate gets converted into radiation, the 
fluctuations in A translate into the fluctuations in the final reheating temperature. 
Note that in our scenario different regions go through slightly different histories, 
as the couplings are different. Thermal equilibrium ensures that we get adiabatic 
fluctuations. 

The key point of our approach is to notice that A is not a constant, and can 
fluctuate in space. The reason why A could fluctuate is that in supersymmetric 
theories as well as in the theories inspired by superstrings, the effective couplings 
are not constants, but rather functions of the scalar fields in the theory. In the early 
universe, these scalar fields can assume different values in different regions. We shall 
come back to this issue below. 

Assuming that A is a stochastic variable, the new mechanism works as follows. 
The decay rate of the inflaton field is 



where m is the inflaton mass at its minimum. We shall assume that V is less than the 
Hubble parameter during inflation. If this were not the case and V » H then the 
reaheating would be instantaneous and all the energy stored in the inflaton would 
immediately be converted into radiation. In that case the reheating temperature 
would be independent of Y. 

When r < H during inflation, the reheating temperature is roughly 



\<pqq 



(1) 



T ~ A 2 m 



(2) 
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If A fluctuates, the corresponding fluctuations in the reheating temperature are: 



$ Tr ST 5\ 

~T R f A" (4) 

The fluctuations are completely determined by the fluctuations of V not of <fi. 

In order to understand in more detail how our mechanism works we will consider 
a toy model. We write down the full set of perturbation equations in a later section. 
Let us follow the cosmological evolution in two different domains of the Universe in 
which A takes two different values equal to Ai and A2 respectively. For definiteness 
we shall assume that Ai > A2. According to our assumption, right after inflation the 
energy densities in the two domains are equal, and evolve as non-relativistic matter, 

tlmatter = ^matter ~ — o (5) 

or 

where a is the scale factor. The inflaton decay in each domain takes place when 

r = h ~ (6) 

M PL 1 ' 

where if is a Hubble parameter. In the approximation of an instant decay (which 
we adopt for the moment) the reheating temperature is given by (jHJ). Since Ai > A 2 , 
the condition (jUJ) gets first satisfied in the Ai-domain. Thus when 

r\ = Xlm ~ H (7) 

the energy in this domain gets converted into radiation: 

eirad ~ T R1 ~ \\m 2 Mp L (8) 

which subsequently scales as 1/a . Right after this moment the energy densities in 
two domains are still equal. However, in the Ai-domain it is stored in radiation, 
whereas in the A2-domain it is still stored in inflaton oscillations and evolves as 
non-relativistic matter. From this moment on up until reheating in the A2-domain, 
the energy densities in the two domains evolve differently. 
The A2-domain gets reheated when 

T 2 = \\m ~ H (9) 

and gets filled with radiation of an energy density 

e 2 rad ~ ~ \\m 2 M 2 PL (10) 

For equal values of scale factors, the energy in the first domain is redshifted to 

/ \ \ 4 /3 

eirad ~ ( y) X A 2 m 2 M 2 PL (11) 
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The resulting radiation energy densities in these domains are related as 



A \ 4/3 

t~ e 2rad (12) 




The density perturbations are 

5e ra d 5\ 5T 

Thus, during the reheating process, the fluctuations in A are translated in den- 
sity perturbations. The precise calculation of the above model (§4 and Appendix) 
actually gives Se rad /e rad « 0.15T/T. 



2 The origin of the fluctuations in the decay rate 

Let us now discuss why the fluctuations in T could be generated in the first place. 
The reason why coupling "constants" can fluctuate, is that in unified theories, such 
as supersymmetric theories, the couplings are not constant but are determined by 
expectation values of fields in the theory. Many of these fields have very flat poten- 
tials and thus can strongly fluctuate during inflation. 

In general we can think of A (S) as a function of a scalar field S. Let us assume 
that q-s are the ordinary fermions, and their supersymmetric partners. The coupling 
can be expanded in series of S 

\(S) = A (l + ^ + ...) (14) 

Where M is some scale, which may be of order M PL or smaller. Since we want to keep 
our discussion maximally general, we should not specify the value of M. Depending 
on the gauge symmetries of the theory the first (and some of the higher terms) of 
expansion may or may not be absent. We shall assume that S is a light field, with 
mass smaller than the Hubble parameter during inflation. Then S fluctuates during 
inflation, 8S ~ H and these fluctuations translate into fluctuations of A, which after 
reheating translate into the density fluctuations. 

We can define / to be the fraction of the coupling controlled by the fluctuating 
field. If the VEV of S during inflation is (S) then / = (S)/M. In terms of / we can 
write 

T = f Js) (15) 

The observed level of density perturbation implies that ST/T ~ 10~ 5 . The fluctu- 
ations of S are of order, 5S ~ H, so the observed level of fluctuations can be the 
result of (S) » H or of / << I. 
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The fluctuations of S will be Gaussian if (S) » H and non-Gaussian in the 
opposite limit. Thus depending on the VEV of S during inflation we can get fluc- 
tuations with a varying level of non-Gaussianities. Our mechanism for generating 
non-Gaussianities is different from the multifield inflation models (eg. [2JE] and ref- 
erences therein). 

2.1 Example 

Let us consider a simple illustrative example: Minimal Supersymmetric Standard 
Model (MSSM) plus an inflaton 0. The chiral superfield content of the model 
consists of two Higgs doublets (we shall call them h, h) , quark and lepton superfields 
plus their anti-particle superfields (we shall denote them as q, q c ), and an inflaton 
superfield ((f)). Following the usual practice, we shall assume that the inflaton is a 
singlet under the MSSM gauge group. The specific form of the inflaton potential is 
unimportant. All we assume is that it gives us a sufficient amount of inflation, for 
solving the flatness and horizon problems. We shall also assume that the density 
fluctuations created by the inflaton fluctuations are very small. 

At the end of inflation, has to decay and reheat the Universe. Let us ask, 
how could decay into quarks and leptons. Since is a gauge singlet, it has three 
options for transmitting its energy into ordinary particles. 

1) can directly decay into the ordinary particles through the renormalizable 
interaction in the superpotential 

X (f)hh (16) 

In fact due to the gauge symmetry, the only candidates for such a decay are the 
Higgs doublets. 

2) can also decay through non-renormalizable interactions in the superpotential 

+<P^q c qc +<P-^qq c + ■■■ (17) 

where M is some mass scale. This scale can be ~ Mrl or may be a lower scale 
coming from integrating out some heavy particles below Mpl- The precise origin 
of the couping is not important. 1 Analogous, couplings may be included in Kahler 
potential. 

3) Finally, can decay into some heavy exotic particles with MSSM gauge 
charges, which quickly re-scatter and thermalize with the ordinary particles. In 
particular, the effective non-renormalizable couplings of the form (jl7j) can result 
from integrating out such states. 

Let us focus on the second possibility. If during inflation and reheating one of 
the sfermions participating in the the coupling gets a VEV, then through (fTTj) can 

^ome of the above interactions break baryon number, however they cannot lead to proton 
decay if (f> has a zero VEV today. Breaking of baryon number conservation could be even welcome 
for generating baryon asymmetry during reheating. This issue will not be discussed here. 
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experience a direct two-body decay into the g-quanta. The decay rate is regulated 
by an effective couping 

* - § 

Where (S) is the VEV of a scalar component of one of the g-superfields. The 
resulting partial density perturbations coming from the above channel will be 

5e through s ~ ^ (j^j m 2 M 2 PL , (19) 

whereas the total density will be sum over all the energies created in different chan- 
nels, including the direct decays through non-fluctuating couplings (JEJ) 

e d irect ~ X 4 m 2 M 2 PL (20) 

(we ignore the delay in thermalization due to re-scattering of the different MSSM 
species.) Thus, the density perturbations will be given by 



0£ r ad ^through S 



^rad ^through S ^direct 



(21) 



The properties of these perturbations, such as Gaussianity will be determined 
by the balance between the different channels, and on the value of 5S. Let us briefly 
discuss this issue. 

In MSSM S corresponds to one of the many flat directions of the potential. 
These flat directions are flat only in the unbroken SUSY limit, but are lifted by 
SUSY breaking effects, and acquire curvature ~ TeV. 

During inflation (or reheating) the inflationary energy density breaks supersym- 
metry and lifts the flat directions even more, provided H > TeV. So in general, the 
flat direction fields acquire masses ~ H and may acquire big VEVs as well (see |3]). 

The physical reason behind the ~ H 2 curvature of flat directions during inflation 
is very transparent. The point is that a positive vacuum energy density, that drives 
inflation, breaks supersymmetry spontaneously. This breaking is either F or a D- 
type, meaning that the vacuum energy density is dominated by either the F or the 
D terms. In case of, F-breaking, the breaking is universally transmitted to all the 
fields through the couplings in the Kahler potential, e.g., of the following form 

d* x d A 8^S*S (22) 

where is the superfield with a non-zero F-term, and S is a flat direction field. 
Such couplings will be generated both by supergravity (in which case M ~ Mp^) 
and by gauge interactions ^] (in which case M <C Mpi). For instance, in the case 
of a minimal Kahler, the gravity-mediated curvature of the flat directions is 

2 2 \ W ? 

m gravity — 3 + 4 (23) 

m PL 
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777^ 

gauge 



where W is a superpotential. However, in the MSSM model the flat directions are 
usually lifted also by gauge-mediated contribution. For instance, if the inflaton <f> can 
couple to some gauge-charged fields (which is usually the case), the flat directions 
receive a two-loop gauge contribution to the masses [3] 

" (=)'(#)' 

One-loop contributions are also possible. The over-all sign of the contribution and 
their balance is very model dependent. For instance, in vD-term inflation 5J gravity- 
mediated curvature is <C if 2 , and the gauge-mediated contribution is dominant. 

The bottomline of our discussion is that, generically, flat directions are lifted 
during inflation by gravitational and gauge effects. As a result their VEVs are 
shifted. The location of the minima, and effective masses in those minima are very 
model dependent and it is impossible to draw a general conclusion. On the other 
hand construction of a particular model which would accommodate our needs is 
trivial. The condition under which the new mechanism of the density perturbations 
will dominate is that the mass of S during inflation be slightly below H (say ~ 0.1H 
or so). Then H behaves as practically massless during inflation and we have 5S ~ H. 

Corrections of the form (|23J1 would inevitably ruin inflation, and must be avoided. 
Inflationary scenarios that universally suppress the gravity-mediated curvature (|23|) 
of flat directions, e.g., such as the D-term inflation [S], are therefore highly motivated, 
as they protect the flatness of the inflaton potential from dangerous gravity-mediated 
corrections ~ H 2 . In such scenarios our mechanism will generically be operative and 
contribute to the density perturbations. 

Then depending on which is the dominant channel of the inflaton decay, we get 
different pictures in perturbation spectrum. If couplings like (fTTj) are the dominant 
source of reheating, then the density perturbations are given by 

e (S) 

In which case to get the correct level of density perturbations we have to demand 
that (S) ~ 10 5 H, and resulting perturbations will be Gaussian, as in the case of 
an ordinary inflation. 

Another possibility is when the channel dominates. In this case 

5e _6S_ _H_ 

7 ~ M\l ~ MAg ( } 

Since now S can stay near its inflationary minimum, the perturbations can be non- 
Gaussian. 



3 Consequences for model building 

We have argued that our new mechanism for generating density perturbations cre- 
ates fluctuations of order H/M where M is a scale that could be significantly lower 
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than the Plank scale (the origin of M depends on the particular model, it could 
be related to the VEV of a field or with physics at scales intermediate between 
Hubble during inflation and the Planck Scale). In other words our mechanism can 
create the observed density perturbations with a much lower scale of inflation than 
the standard scenario where density perturbations are given by H/MpLt/e, with 
e=(M PL V'/V) 2 . 

If our mechanism is responsible for the observed perturbations then we predict 
a much lower background of gravitational waves than in standard scenarios. The 
amplitude of the gravity wave background is still given by H/Mpi and thus it is 
greatly suppressed relative to density perturbations, of order H/M. 

The relation between the slope of the power spectrum of density perturbations 
and the characteristics of the inflaton potential is different in our scenario, 

dlnH 2 

n — 1 = — our scenario 

a ma 

d\nH 2 /t 

n — 1 = — standard case (27) 

am a 

As a result, both the energy scale of inflation and the characteristics of the in- 
flaton potential would be misinterpreted if our mechanism was operating but the 
observations were analyzed using the standard assumptions. 

Moreover, there is a significant chance that the fluctuations generated by our 
mechanism could be at least slightly non-Gaussian. The fields that determine the 
couplings are not slow rolling during inflation and do not dominate the energy den- 
sity. As a result, their potentials need not satisfy the stringent slow- roll conditions 
the inflaton needs to satisfy. These conditions imply that non-Gaussianities gener- 
ated in standard one-field models of inflation are very small, probably unobservable 
with current techniques The fields determining the couplings could have for 
example larger self couplings (such would be, for instance, the quasi-flat directions 
stabilized by quartic Yukawa couplings in the superpotential) leading to observable 
non-Gaussianities. 

The slow roll condition on the kth. derivative of the inflaton potential (V^) is, 

^M^VV < L (2g) 

For k = 3 for example, this leads to the constraint that / H < H / Mpie 1 ^ 2 ~ 
10 -5 . In our scenario couplings much larger that this are allowed. The cubic coupling 
over H is roughly given by m 2 s / {S)H ~ H/ (S) so if (S) is not much larger than H 
significant non-Gaussianities are possible. 

The bound on the three point function from WMAP can roughly be translated 
into V™/H < 10~ 3 [Zj. On the other hand to match the correct level of density 
perturbations f5S/ (S) ~ 1CT 4 , (as we will show in the next section the solution 
of the full equations imply that the potential fluctuations are of order 1/9 SF/F so 
we set the combination fSS/ (S) to 10 -4 rather than 10~ 5 ) . This means that when 
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/ ~ 0.1 the non-Gaussianities saturate the WMAP bound. A more detailed analysis 
of the constraints set by WMAP is left for a future publication. 

4 Perturbation Equations 

In this section we will explicitly write the relevant perturbation equations in the 
conformal gauge [3IH]- We will consider two fluids, the inflaton field which for 
simplicity we will treat as non-relativistic matter and radiation. 

First we write the equations for the background evolution. The background 
metric is flat FRW universe, ds 2 = a 2 (dr] 2 — dx 2 ), rj is the conformal time, a' /a = 7i 
is the conformal Hubble constant, with ' = d/drj. The unperturbed Friedmann 
equations are 

e' m = -3He m - Tae m , (29) 
e' r = -me r + Tae m , (30) 
H 2 = (8/3)nGa 2 (e r + e m ). (31) 

These describe a matter-dominated universe at small time and a radiation- 
dominated universe at large time. We solved them assuming that expansion rate was 
if* and x = e m /(e m + e r ) = 1 at the beginning. The fraction x remains around one 
until it suddenly drops to x = at the time when H ~ T. Before that time H 2 oc a~ 3 
and after that H 2 oc a~ A . Here H = 7i/a is the standard Hubble parameter. 

We now consider super-horizon perturbations induced by perturbations in the 
decay rate V. We use the conformal gauge with perturbed metric ds 2 = a 2 ((l + 
2&)di] 2 — (1 — 2Q)dx 2 ). The perturbed Einstein and matter equations describing 
forced super-horizon perturbations are 

H& + H 2 § = -(4/3)nGa 2 (e m 6 m + e r 6 r ), (32) 

5' m = 3$'-ra(5 r + $), (33) 

S' r = 4$' + (e m /e r )Ta(5 r + $ + 5 m - 5 r ). (34) 

Here 5 m = 5e m /e m , 5 r = 5e r /e r , 5 r = 5T/T. 

We solved equations (|32H34|) numerically. The solution for the gravitational 
potential is shown in figure ^ The gravitational potential directly gives the level of 
Microwave Background anisotropics which are simply proportional to $ (see jH] for 
details on the calculation of CMB anisotropies in the conformal gauge). As we have 
estimated before, the potential fluctuations are of order <5r/r. More specifically $ m 
1/9 ST/T with only a weak dependance of T /H*, provided T/H* is relatively small. 
In the limt T/H* 0, one obtains (Appendix) the exact result $ = (1/9)ST/T. 
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Figure 1: Evolution of the gravitational potential in units of 5T/T as a function of 
iV = In a. The curves correspond to if*/T = 5, 10,20 (solid, dotted, dashed). 

5 Conclusions 

We have proposed a new way in which density perturbations could be created during 
inflation. The density fluctuations produced by the new mechanism could easily 
dominate over the ones produced in the standard scenario. Our mechanism allows 
for the energy scale of inflation to be significantly lower than usually assumed. 
The fluctuations produced are also adiabatic and nearly scale invariant, although 
the relation of the spectral slope to the inflation potential is different from the 
usual. If the scale of inflation is lower, we generically predict a lower background of 
gravitational waves. 

Our mechanism has no trouble generating non-Gaussianities that could be ob- 
served by future experiments. If such signal were to be found it would give additional 
clues as to the details of the inflationary model and how it fits in the general frame- 
work of particle physics. 
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A Appendix: Isomorphism between free and in- 
duced perturbations. Exact solution. 

Equations for free perturbations are obtained just by dropping the terms with 5r 
from equations for induced perturbations 



H& + H 2 $ = -(A/3)nGa 2 (e m 5 m + e r 5 r ), (35) 

S' m = 3$' - Ta$, (36) 
8' r = A& + (e m /e r )Va(^ + 5 m -6 r ). (37) 

Let S m , 5 r be a solution of the free perturbation equations (|35H37|) . Then 
$i = $ — Sr, S m i = 5 m + 25r, 5 r \ = 5 r + 2<5 r is a solution of the forced perturbation 
problem. We are interested in the purely forced mode, when $ x = 5 ml = 5 r \ = at 
7] = 0. The purely forced mode can be obtained from the free solution which satisfies 
$ = 5r, 8 m = —25r, S r = —25r at rj = 0. This is just an adiabatic perturbation 
in the matter-dominated universe. That 5 r ^ at rj = is irrelevant, because the 
radiation energy is small at small rj. 

We will show below that the free adiabatic (not growing for 77 — ^ 0) perturbation 
satisfies the well known matching condition [TU] 

$(00) = = 10 (3g) 

$(0) l+p r 9' 1 ' 

Here the indices p characterize the expansion rate of the universe, a oc t p . For matter 
domination p m = 2/3, and for radiation domination p r = 1/2. 

The problem of calculating the forced metric perturbation is thus solved. The 
fluctuating reaction rate 5T induces a super-horizon potential perturbation 

(30) 



It remains to prove (|38j). The free perturbation equations (|35H37|) can be written 

as 

(3e m + 4e r ) 2 C = e m S, (40) 

s> = - ^ + ^ r ns, (4i) 

where ( (the Bardeen parameter) and S are 

C = $ - ^±^ ; (42) 
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S = 7ie r {2>8 r — 45 m ) + Ta(e m 5 m + e r 5 r ). (43) 

There are two scalar modes of free perturbations. The mode S 7^ is obtained from 
(141)1 . The other mode is S = 0, ( = const. One can show that the latter mode 
describes adiabatic perturbations at 77 — > 0. We can therefore use ( = const to prove 
the matching condition 
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